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INTRODUCTION 

Hyper-accretion disks form when a neutron star merges with another compact object, 
neutron star or a black hole. Recent numerical simulations of mergers 01 0, 1] are fully 
relativistic and show how most of the mass of the binary system disappears behind the 
event horizon in about 10 ms, leaving a rotating debris disk around the black hole. The 
mass of this centrifugally supported disk is m ~ 0.01 — 0.1 M , similar to what was 
found in previous non-relativistic simulations (e.g. [0, H]). The ensuing disk accretion 
is not followed by the merger simulations. It is established on a viscous timescale 
^visc ~ 0. 1 s as the inner parts of the disk relax to a quasi-steady state; this relaxation was 
studied numerically in [a [ZD. Most but not all of the disk is accreted on the timescale 
/vise* releasing an energy comparable to mc 2 and emitting copious neutrinos. This disk 
has an accretion rate M ~ (M /s) (m/0.1M Q ) (j visc /0.1 s) _1 . 

Similar neutrino-emitting disks may form during the core-collapse of massive stars 
if the stellar material has a sufficient angular momentum [[8|, §]. These hypothetical 
objects are often called "collapsars." After the formation of a central black hole of mass 
M ~ 3 M collapsars develop an accretion disk that is fed by the continually infalling 
stellar material. The high accretion rate M ~ 0. 1 M s _1 is sustained for ~ 10 s (the core- 
collapse timescale). Recent relativistic MHD simulations of this accretion show how the 
black hole could accumulate a large magnetic flux and create jets via the Blandford- 
Znajek process Qp|l . 

The studies of hyper-accretion disks are greatly stimulated by observations of cos- 
mological gamma-ray bursts (GRBs, see lijj] for a review). Hyper- accretion is ex- 
pected to produce hyper-jets on a timescale ~ 0.1 — 10 s. If a fraction £j et of the ac- 
cretion power Mc 2 is channeled to a relativistic jet, it leads to an explosion with energy 
£j et ~2x lO 51 (M acc /M )(£j et /lO~ 3 ) erg, where M acc is the mass accreted through the 
disk. The energy and duration of the jet is consistent with GRB observations. 

Hyper- accretion disks are markedly different from normal accretion disks in X-ray 



binaries and AGN. Their optical depth to photon scattering is enormous and radiation 
is trapped inside the disk, being advected by the matter into the black hole. However, 
the disk can be efficiently cooled by neutrino emission. Significant neutrino losses can 
occur when M > 10~ 3 M s _1 and make the disk relatively thin and neutron rich. 

The accreting black hole is expected to have a significant angular momentum, because 
it forms from rotating matter and is further spun up by accretion. The black-hole spin 
helps the jet formation through the Blandford-Znajek process. It also affects the black- 
hole spacetime in such a way that the disk extends to smaller radii and the overall 
efficiency of accretion significantly increases. For example, the inner radius of the disk 
around a maximally -rotating black hole (spin parameter a = 1) is reduced by a factor of 
6 compared with the Schwarzschild case a = 0. This leads to a higher temperature and a 
higher neutrino intensity above the disk, increasing the rate of neutrino annihilation into 
e ± pairs. Therefore, disks around rapidly spinning black holes can deposit an interesting 
fraction of their energy into the e ± plasma outside the disk and facilitate the formation 
of ultra-relativistic jets. 

The size of a hyper-accretion disk depends on the specific angular momentum of the 
accreting matter, /, which is modest in neutron-star mergers and probably even smaller in 
collapsars. The accretion flows in collapsars are quasi- spherical and may form a special 
"mini-disk" that is not supported centrifugally and instead accretes on the free-fall time. 
Larger / leads to standard viscous accretion, which leaves a relict disk carrying the initial 
angular momentum of the accreted matter. The relict disk gradually spreads to larger 
radii, and its late evolution may be relevant to the post-burst activity of GRBs. 

VISCOUS DISKS 

As matter spirals into the black hole, it is viscously heated: the gravitational energy 
is converted to heat. The heat is distributed between nuclear matter, radiation, and e ± 
pairs, in perfect thermodynamic equilibrium. In particular, the equilibrium e ± population 
is maintained. As discussed below, electrons are mildly degenerate in neutrino-cooled 
disks, which affects the density. The equilibrium microphysics is determined by only 
three parameters: temperature T, baryon mass density p, and electron fraction Y e (equal 
to the charged nucleon fraction). Other parameters — e.g. the electron chemical potential 
jl e and density of pairs, n± — are derived from T, p and Y e . At radii r < 10 8 cm, 
temperature and density are high enough to maintain the nuclear statistical equilibrium, 
which determines the abundances of all nuclei. Nuclear matter in the disk is dissociated 
into free nucleons n and p inside radius r a = (40— I00)r g where r g = 2GM/c 2 . The 
temperature at this radius is kT < 1 MeV. 

Neutrino cooling is significant in the inner region where kT > 1 MeV. The far domi- 
nant mechanism of neutrino emission is the capture onto nucleons: 

e~ +p^n + v e , e + +n -> p + V e . (1) 

The escaping neutrinos not only cool down the disk — they also change its electron 
fraction Y e if the emission rates of v and v are not exactly equal. The first goal of the 
disk modeling is to find T, p, Y e , and self-consistently evaluate the neutrino losses. 



Modeling disk accretion 



Accretion is quasi-steady in the inner region of the disk where t w [ sc is shorter than the 
timescale of M evolution. Viscosity in accretion disks is caused by MHD turbulence, 
sustained by the magneto-rotational instability lfl2ll . It creates an effective kinematic 
viscosity coefficient v, which may be related to the half-thickness of the disk, H, and 
sound speed, c s : V = (Xc s H, where a ~ 0.01 —0.1 is a dimensionless parameter. A 
number of works studied hyper- accretion disks with this traditional parameterization of 
viscosity lfl3l [Til [l5l [l6l 17 . 18 . 191. Two works constructed models in full general 



relativity and studied disks around spinning black holes II13L Il7|l . The following list 
highlights advances and limitations of the current models : 

- The model is fully relativistic: the disk dynamics is calculated in Kerr metric. 

- Hydrodynamic equations are solved with the vertically-integrated a prescription. The 
equations include radial transport of heat and lepton number. 

- Local microphysics is calculated exactly: nuclear composition, electron degeneracy, 
neutrino emissivity and opacity etc., using the equilibrium distribution functions for 
all species except neutrinos. Neutrinos are modeled separately in the opaque and 
transparent zones of the disk, matching at the transition between the two zones. 

- The model provides only vertically averaged T, p, and Y e (which are approximately 
equal to their values in the midplane of the disk). The diffusion of neutrinos in 
the opaque zone is treated in the simplest one-zone approximation (using escape 
probability). This gives a good approximation for the energy losses, however, does 
not give the exact neutrino spectrum emerging from the opaque zone. 

The vertically-integrated approximation provides no information about the vertical 
structure of the disk and its corona. The vertical structure may be eventually understood 
with global 3D time-dependent MHD simulations that include energy losses, although 
the results of such simulations generally depend on the assumed initial magnetic con- 
figuration. The behavior of magnetic field on large scales is coupled to the local turbu- 
lence cascade that extends to scales <C H. The microscopic magnetic Prandtl number for 
hyper-accretion disks has been recently estimated in Q20Q. The models discussed here are 
aware of the MHD issues only through the value of a. They are computationally much 
cheaper than MHD simulations and allow one to study the disk in a broad range of M 
and a. 

Vertically-integrated disks are described by ID equations that express conservation of 
baryon number, energy, and momentum (angular and radial) in Kerr spacetime (see [21] 



for a review). The full set of these equations can be solved |22j], and the solutions show 
that the deviation from circular Keplerian rotation is small (< 10%) even when the disk 
is strongly advective (i.e. when the released heat is transported radially without losses)Q 
Thus, the angular velocity of the disk can be approximated by its Keplerian value Qk- 
A small radial velocity is superimposed on this rotation: u' = —aS l c s (H/r), where 



A strong reduction of £2 below £2k can occur in the limit of a large steady disk with no cooling. This 
limit does not apply to hyper-accretion disks which are transient and have a moderate radius. 
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FIGURE 1. Left panel: Contours of the equilibrium Y e (T, p ) on the T-p plane for v-transparent matter. 

The electrons become degenerate near the dashed line given by kT deg = hc(p/m p ) 1 ^ = 7.7 'p^ 3 MeV. 
The Y e contours are calculated assuming that the nuclear matter is dissociated into free nucleons; they 
are invalid in the shaded region where matter is dominated by composite nuclei. The "neutronization 
line" Y e = 0.5 is given by kT n = 33p ll MeV. Right panel: The equilibrium Y e (T,p) for v-opaque matter 
with neutrino chemical potential ji v = 0. The free-nucleon region is the same as in the left panel. The 
calculation of the equilibrium Y e is now extended into the region of composite nuclei. The neutronization 

line Y e = 0.5 is given by kT n = 23.\p{( 2 MeV. 



c s = (P/p) 1 / 2 is the isothermal sound speed, H is the half-thickness of the disk; S(r) 
is a numerical factor determined by the inner boundary condition [fivh . This description 
of the velocity field in the disk is a good approximation everywhere except in the very 
vicinity of the inner boundary where \u r \ exceeds c s . 

In contrast to accretion disks in X-ray binaries and AGN, there is one more conserva- 
tion law that must be taken into account: conservation of lepton number, 
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Here A^ v and /Vy are the number fluxes of neutrinos and anti-neutrinos per unit area (from 
one face of the disk), n v and ny are the number densities of neutrinos and anti-neutrinos 
inside the disk. This equation determines Y e , which is related to the neutron-to-proton 
ratio by Y e = (n n /n p + 1) _1 and greatly affects the rate of neutrino cooling. 

In the models solved in [fivh and shown below, Y e is calculated using Eq. ©. Note, 
however, that throughout most of the neutrino-cooled disk, the right side of Eq. © is 
small compared with each of the two terms on the left side, and Y e is nearly equal to the 
local equilibrium value such that N V ^N V . This equilibrium Y e is determined by the local 
temperature and density and found for both neutrino-opaque and neutrino-transparent 
matter EBB]- It is shown in Fig. 1. 
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FIGURE 2. Boundaries of different regions on the r-M plane for disks around a black hole of mass 
M = 3 M and spin parameter a = 0.95. Neutrino cooling is inefficient in the shaded region below the 
"v-cooled" curve and above the "trapped" curve. The shaded region marked "unstable" is excluded: the 
steady model is inconsistent in this region because of the gravitational instability. The disk extends down 
to the marginally stable orbit of radius r ms w r g where r g = 2GM/c 2 . Left panel: Disks with viscosity 
parameter a — 0.1. Right panel: Disks with viscosity parameter a = 0.01. (From (lTn .) 



Overview of disk properties 

Hyper- accretion disks have several zones separated by the following characteristic 
radii: 

1. Radius r a where 50% of a-particles are decomposed into free nucleons. The 
destruction of a-particles consumes 7 MeV per nucleon, which makes the disk 
thinner. 

2. "Ignition" radius r lgri where neutrino emission switches on. At this radius, the mean 
electron energy becomes comparable to (m„ — m p )c 2 , enabling the capture reaction 
e~ + p — > n + v. Then neutrino cooling due to reactions ([I]) becomes significant, 
further reducing the disk thickness H jr. 

3. Radius r v where the disk becomes opaque to neutrinos and they relax to a thermal 
distribution. The disk is still almost transparent to anti-neutrinos at this radius. 

4. Radius where the disk becomes opaque to anti-neutrinos, so that both v and V 
are now in thermal equilibrium with the matter. The disk is still cooled efficiently 
at this radius since v and v diffuse and escape the flow faster than it accretes into 
the black hole. 

5. Radius r tl - where neutrino diffusion out of the disk becomes slower than accretion, 
and neutrinos get trapped and advected into the black hole. 

The different zones of the disk are shown on the M — r diagram in Fig. 2. In addition, 
this figure shows the zone of gravitational instability. 



Three characteristic accretion rates can be defined: Mi gn above which the disk is 
neutrino-cooled in the inner region, M paque above which the disk is opaque to neutrinos 
in the inner region, and M ttap above which the trapping of neutrinos occurs in the inner 
region. The dependence of M lgn , M paque, and A/ tra p on a is well approximated by the 
following power laws DJ/ZD, 
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The normalization factors K depend on the black hole spin a. For a = 0.95 they are 
K iga = 0.021 M^s" 1 ,^ 



; ^ ^opaque = 0.06 M Q s 1 , _K t rap = 1 -8 M s 1 and for a = they are 
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To complete this short guide to quasi-steady viscous disks, Figs. 3-4 show T, p, Y e , 
and H/r for a model with M = 0.2 M s -1 , M = 3 M , and a = 0.95. 
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FIGURE 3. Disk with M = 0.2 M s~ 1 around a black hole of mass M = 3 M Q and spin a = 0.95. Three 
models are shown with viscosity a = 0.1, 0.03, and 0.01. Radius is measured in units of r g = 2GM/c 2 = 
10 km. Left panel: Temperature in units of m e c 2 . Right panel: Mass density. (From {vh .) 
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FIGURE 4. Electron fraction Y e and thickness of the disk H/r for the same three models as in Fig. 3. 
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FIGURE 5. Contributions to total pressure P from baryons, Pj,, electrons and positrons P e = P e - + P e + , 
radiation Py, and neutrinos P v + P v for the accretion disk with M = 0.2 M Q s around a black hole of 
mass M = 3 M Q and spin a — 0.95. Left panel: Model with viscosity parameter a 
Model with a = 0.01. (From 0.) 
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The main properties of the neutrino-cooled disk (i.e. at r < ri gn ) may be summarized 
as follows S. 

o The disk is relatively thin, H /r ~ 0. 1 — 0.3, especially in the inner region where most 
of accretion energy is released. 

o The v-cooled disk is locally very close to /3 -equilibrium, N v « Ny. In particular, the 
relation between p, T, and Y e calculated under the equilibrium assumption (Fig. 1) is 
satisfied with a high accuracy. 

o Degeneracy of electrons in the disk significantly suppresses the positron density n e +. 
However, the strong degeneracy limit is not applicable — the disk regulates itself 
to a mildly degenerate state with ll e /kT = 1—3. The reason of this regulation is 
the negative feedback of degeneracy on the cooling rate: higher degeneracy }l e /kT 
— > fewer electrons (lower Y e ) and positrons (n e +/n e - ~ A r ) weaker neutrino 
emission — > lower cooling rate — > higher temperature — > lower degeneracy. 

o Pressure in v-cooled disks is dominated by baryons, P ~ P/, — (p/m p )kT, most of 
which are neutrons (Fig. 5 shows contributions to pressure for two sample models). 

o All v-cooled disks are very neutron rich in the inner region, with Y e ~ 0. 1 or lower. 



NEUTRINO ANNIHILATION AROUND THE DISK 

The emitted neutrinos and anti-neutrinos can collide and convert to e ± , thereby deposit- 
ing energy feol . The emission of tau and muon neutrinos is negligible II 1711 . so only 
reaction v e + v e — > e~ + e + can be considered. Its cross section (assuming center-of- 
momentum energy 3> m e c 2 ) is given by a vv ~ 3.3 x 10 _45 (pv-pv) 2 (PvP9) _lcm2 wnere 
p v , Pv> are the 4-momenta of V e and v e , expressed in units of m e c. The cross section is 



small and only a small fraction e ~ 10 3 — 10 2 of the total neutrino luminosity L con- 
verts to e ± plasma. Nevertheless, this energy may be sufficient to drive a relativistic jet 
(or help the formation of a magnetically-dominated jet) since it occurs above the disk 
where the mass density is relatively low, especially near the rotation axis ||27|l . 

Neutrinos emitted by the disk follow null geodesies in Kerr spacetime. The efficiency 
e of their annihilation can be calculated numerically by tracing the geodesies, evaluating 
the local energy deposition rate q V y [erg s _1 cm -3 ] everywhere around the black hole, 
and then integrating over volume to obtain the net energy deposition rate E vv (energy at 
infinity per unit time at infinity). The neutrino emission and annihilation is concentrated 
near the black hole, where accretion is expected to be quasi-steady. E V y depends on four 
parameters that specify the steady disk model: M, a, M, and a. 

The energy deposition rate E V y was estimated in |fl3|l . approximating geodesies by 
straight lines. Fully relativistic calculations were made for several toy models, in par- 
ticular for disks or tori with uniform temperature or other arbitrary distribution of tem- 
perature or entropy (see [|280 and refs. therein). Recently, the relativistic calculation for 
a realistic disk around a spinning black hole has been done (Zalamea & Beloborodov, 
in preparation). q V y has been obtained everywhere around the black hole, including its 
ergosphere, and the dependence of E V y on M, a, M, and a has been determined. 

Besides tracing the geodesies, this calculation involves a model for neutrino and anti- 
neutrino spectra emitted by the disk. Using the results of U7\\ . it is straightforward to 
evaluate the spectra from the transparent zone of the disk. It is more difficult to find the 
spectrum that emerges from the opaque zone, because the neutrino transport in this zone 
depends on the unknown vertical distribution of viscous heating. Various assumptions 
may be made about this distribution 000,10, including a strong heating 
of the magnetic corona above the disk. Note that the corona of a hyper-accretion disk 
is always in thermodynamic equilibrium, and its temperature T c is determined by its 
thermal energy density U c = Uy + U e ~ 3a r T^ where a r = 7.56 x 10~ 15 erg cm~ 3 K -4 
is the radiation constant. U c is generally smaller than the energy density inside the disk. 
Therefore, relocating the heating from the disk to its corona cannot significantly increase 
the energies of emitted neutrinos. 

Fortunately, a robust estimate can be obtained for E V y in spite of the uncertainty in 
the vertical structure of the disk. It is easy to see that all detailed models of neutrino 
spectrum formation must predict practically the same rate of vv annihilation above a 
neutrino-cooled disk. For such a disk, neutrinos carry away a fixed energy flux F m F + 
where F + ~ 3MQ.^S(r) /Sn is the rate of viscous heating. Therefore, transfer models 
that predict a higher average energy of emitted neutrinos, E w , must also predict a 
lower number density of the neutrinos above the disk, n ~ F /E av c E~ l . Since the 
annihilation cross section o V v 01 E% y , one finds that the reaction rate h V y ~ c(7 V v«v«v 
is independent of E ay . The energy deposition rate q V y ~ n V yE ay is proportional to E aw . 
It cannot be changed without a substantial change in temperature (or electron chemical 
potential) of the neutrino source, which would require a huge change in energy density 
and therefore is hardly possible. 

Detailed calculations confirm that £Vv weakly depends on the details of the vertical 
structure and v, V transfer in the disk. Consider two extreme models for the opaque zone. 
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FIGURE 6. Total energy deposition rate due to vv annihilation outside the black-hole horizon, £Vv? as 
a function of the disk accretion rate, M. The two characteristic accretion rates M; gn and A/trap depend on 
the viscosity parameter a — 0.1 (see eq.[3); a = 0.1 is chosen in this figure. The black hole is assumed to 
have mass M — 3 M . E V y strongly depends on the spin parameter of the black hole; the numerical results 
are shown for two cases: a = (triangles) and a = 0.95 (squares). The uncertainty in the vertical structure 
of the accretion disk leads to a small uncertainty in £Vv as illustrated by two extreme models: Model A 
(open symbols) and Model B (filled symbols), see the text for details. The results of both models are well 
approximated by simple Model C that is given by Eq. (0]i and also shown in the figure, by solid line for 
V-opaque disks (M > M opa q Ue ) and by dashed line for v-transparent disks (M < Mopaque)- (From Zalamea 
& Beloborodov, in preparation.) 



Model A: Neutrinos v e and v e are emitted with the same distributions as found inside the 
disk (same temperature T and chemical potential /i v ). The distribution normalization is, 
however, reduced compared with the thermal level inside the disk, so that the emerging 
emission carries away the known energy fluxes F v and Fy that are found in [|17ll . 
Model B: Neutrinos are emitted with a thermal Fermi-Dirac spectrum with chemical 
potential /i v = and the effective surface temperature T e ff . The effective temperature is 
defined by (7 /8)oT e f f = F , where F = F v +F v is the energy flux from one side of the 
disk and o = a r c/A is Stefan-Boltzmann constant. The factor 7/8 is determined by the 
difference between Plank and Fermi-Dirac distributions and the fact that the disk emits 
only v e and v e — the emission of other neutrino species is negligible. 

1/4 

The temperature T inside an opaque disk exceeds T e ff by factor T /T e ff ~ t v when 
the neutrino optical depth t v 3> 1 . The neutrino chemical potential is modest and, ap- 
proximately, E &w T . Hence, even the extreme Model A gives a moderate enhancement 
of E V v, by the factor T/T e ff. Zalamea & Beloborodov (in preparation) calculated E V y 



in both Models A and B. The results of numerical calculations are shown in Fig. 6 and 
demonstrate that the difference between the two models is indeed small. 

It is instructive then to consider Model C: same as Model B except that F = F + is 
assumed at all radii. The assumption is clearly incorrect outside the region r tr < r < r lgn . 
Nevertheless, this simplest model gives a good approximation to E V y in a broad range 
Mj gn < M < Mtrap (Fig. 6)H Note that E V y in Model C is explicitly independent of Of. 
However, the range of M where Model C is applicable depends on a (Eq.0). 

The scaling of E V y with M is easy to evaluate analytically. The effective surface 
temperature is related to F m F + by F M. The neutrino number density 

above the disk is proportional to 7^ oc M 3 / 4 . The annihilation cross-section o V y 

(assuming kT s ff > m e c^\ Hence the energy deposition rate E V y °^ 7^ M 9 / 4 , and one 
can write 

M ^ 9 / 4 



E vv = Eq(o) [j^— i J > M ign < M < M trap . (4) 

The normalization factor Eq depends on the black hole spin a and must be calculated 
numerically. For example, a = 0.95 gives Eq « 10 52 erg s _1 , which implies the efficiency 
£ = E V y/L w O.O5(M/M s -1 ) 5 / 4 . It is much larger than the corresponding value for a 
non-rotating black hole, by two orders of magnitude. 

The strong dependence of E V y on a may be seen from the following rough estimate. 
The neutrino luminosity L peaks at r pea ]< that is a few times the inner radius of the disk, 
7ms («) — the marginally stable orbit, which is determined by a. The luminosity depends 

on r ms approximately as r m l, and T e f f at r peak scales as {L/r^ e2ik ) l / A rj'' 4 . The energy 

deposition rate q V y scales as T e 9 ff and E V y scales as r^q V y, which yields E V y <=c r ms 15 ^ 4 . 
Then the reduction in r ms by a factor of 3 (as a increases from to 0.95) gives a factor 
of 60 in E V y. This estimate neglects the fact that the gravitational bending of neutrino 
trajectories is stronger for smaller r ms . Stronger bending implies a larger average angle 
between neutrinos and leads to an additional enhancement of E vv . Therefore, a steeper 
dependence of E V y on r ms is found in numerical simulations. A simple power law 
E V v oc r~ 4 7 is an excellent approximation to the numerical results for < a < 0.95 
which corresponds to r g < r ms < 3r g (Zalamea & Beloborodov, in preparation). 

Note that E V y is defined as the total energy deposition rate outside the event horizon 
(including the ergosphere). A significant fraction of the created e ± plasma must fall into 
the black hole, and only the remaining fraction of E V y will add energy to the jet. This 
fraction depends on the plasma dynamics outside the disk, which is affected by magnetic 
fields and is hard to calculate without additional assumptions. 



2 The deviation of Model B from Model C at M ~ 2M; gn is caused by the overshooting of F above F + , 
which happens just inside of rj 2n (see Fig. 7 and 16 in lfl7ll ). As hot matter accretes into the neutrino- 
cooled region, its stored heat is quickly emitted with F /F + reaching ~ 2 at r ~ Hgn/2. For disks with 
Hgn/2 ~ a few r g , this leads to the enhancement of E V y by the factor ~ 2 9 / 4 compared with Model C that 
assumes F = F + . 
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FIGURE 7. The shaded region shows the range of angular momenta /o that form a mini-disk with 
insufficient centrifugal support, leading to accretion on the free-fall timescale. Iq is defined as the angular 
momentum of the accretion flow in the equatorial plane; angular momentum decreases toward the polar 
axis (shells r = const are assumed to have a uniform angular velocity Q. <C £2k at r 3> r g ). (From Zalamea 
& Beloborodov, in preparation). 



LOW-ANGULAR-MOMENTUM DISKS IN COLLAPSARS 

The quasi- spherical accretion flows in collapsars create a centrifugally supported disk 
if the circularization radius of the flow is sufficiently large, r c i rc ~ 10r g . A smaller 
disk may not be centrifugally supported and then will accrete on a free-fall timescale 



Q33D - h accretes so fast (super-sonically) that the effects of viscosity can be neglected. 



A steady model of this "mini-disk" was constructed in [33] and 2D time-dependent 
hydrodynamical simulations were performed in Ii3411 . 

The mini-disk can be thought of as a caustic in the equatorial plane of a rotating 
accretion flow. It absorbs the feeding infall, and this interaction releases energy, making 
the accretion radiatively efficient. With increasing angular momentum, the size of the 
disk grows up to \4r g c, at which point the centrifugal barrier stops accretion, so that it 
can proceed only on a viscous timescale. Thus, the mini-disk model fills the gap between 
two classical regimes of accretion — spherical (Z < r g c) and standard accretion disk 
(/ 3> r g c) — and is qualitatively different from both. 

The calculations of Q33Q were limited to the case of a Schwarzschild black hole. 
Recently, the model has been extended to the case of a Kerr black hole (Zalamea & 
Beloborodov, in preparation). Fig. 7 shows the range of angular momenta that lead 
to mini-disk formation around a black hole of spin < a < 1. The critical angular 
momentum for viscous disk formation sets the maximum radius of a mini-disk. This 
radius is « 14r g = 2SGM/c 2 for a = and « 5r g for a = 0.95. 



The mini-disk is sandwiched by shocks through which the infalling matter enters the 
disk. The model of 113 311 assumes that the cooling timescale of the postshock material is 
sufficiently short, shorter than the accretion timescale in the disk. Let us check if this 
assumption can be valid for collapsars. The postshock mass density and energy density 
can be estimated as 



p ~ 10 9 %M 0A M 3 2 ( y\ 1 gem" 3 , t/~0.2pc 2 iy ) . (5) 



g 



Here £ ~ 2 is the compression in the shock (note that the shock dissipates only the 
normal component of the infall velocity); M .i = M/0.1M© s~ ! andM 3 =M/3M . The 
postshock matter has \l e jkT < 1, and U is dominated by radiation and e ± pairs, which 

implies U ~ 3a r T 4 . The postshock temperature is then T ~ 5.3 x 10 10 p^ 4 (r/r g )~ 1 / 4 K. 
Disintegration of nuclei in the shock consumes only ~ 10~ 2 of the energy released 
at r ~ r g , so the postshock matter can be cooled only by neutrino emission. Neutrino 
emission is dominated by two processes: (i) capture reactions CD provide cooling rate 
q c 9 x 10 32 r ] 6 1 p9 erg cm~ 3 s _1 , and (ii) e ± annihilation e + + e~ — > V + V provides 
q± « 3.6 x lO 33 ^ erg cm~ 3 s _1 (see e.g. n35\\ ). This gives 



* V/g/ cm s \ r gJ cm s 

Approximating the total q = q c + q± ~ g c , one finds 

where ? a cc ~ 10- 4 (r/r g )- 3 / 2 M 3 s. The mini-disk is neutrino-cooled if q t acc /U ^ 1, 
which requires a high accretion rate, comparable to M s _1 . Thus, only high-M mini- 
disks are sandwiched by radiative shocks that stay near the equatorial plane in the 
innermost region of the accretion flow. A large neutrino luminosity, up to ~ O.lMc 2 , 
is produced by such disks. 

For smaller accretion rates, the postshock matter is unable to cool on the free- 
fall timescale, and the neutrino luminosity from the inner region is suppressed by 
the factor qt acc /U < 1. Then a hot low-angular-momentum bubble must grow around 
the black hole. Mass flows into the bubble through the shock front that expands to 
r ^> r g . Such a bubble is observed in low-M simulations in j34 \. It resembles the bub- 
ble around viscous disks in the models of [@, except for a slower rotation, less 
centrifugal support, and faster accretion. The shock expansion can be stopped when 
it approaches ~ 40r g = 4 x 10 7 M3 cm. Up to this radius, the postshock temperature, 

T ~ 7 x 10 9 Mo^ 4 (r/40r g ) -5 / 8 K, is high enough to disintegrate nuclei at the density 
p ~ 4 x 10 6 Mo.l( r / 40r g )~ 3 / 2 g cm~ 3 . As the shock expands to 40r g , its energy de- 
creases to GMm p /r ^ \2(r/A0r g )~ l MeV per nucleon, and a large fraction of this en- 
ergy is consumed by disintegration (8 MeV per nucleon); therefore, the shock stalls. 



SPREADING OF VISCOUS DISKS AND NUCLEAR BURNING 

Formation of a viscous disk with r c i rc > 10r g implies that most of the angular momentum 
of accreting matter will be stored outside the black hole, in a viscously spreading ring. 
When matter supply to the disk stops, accretion will proceed from this ring. At any time 
t, the characteristic size of the ring R(t) is where its mass peaks. Alternatively, R can be 
defined by J = (GMR) l / 2 m, where J is the angular momentum carried by the disk and 
m is its mass. 



Spreading of merger disks 

Immediately after the merger, the characteristic size of the debris disk is Rq ~ 10 7 cm, 
and its initial mass mq may be as large as ~ O.1M | Its viscous evolution starts on a 
timescale to = ( ocQk) ~ 1 (H/r) ~ 2 ^ 0. 1 ( a/0. 1 ) ~ 1 s, with accretion rate M ~ m o Ao that 
can exceed 1 M Q s . The disk is initially hot and v-opaque (cf. Fig. 2); its nuclear 
matter is composed of free nucleons n and p. 

The initial accretion phase lasts ~ to- Following this stage, the disk mass m(t) is 
reduced and its radius R(t) grows to conserve the angular momentum, J = m(GMR) l l 2 
const, which implies m R~ l l 2 . Several important changes occur in the disk as it spreads 
toi?~ 10 2 r g ^ 10 8 cm: 

o Temperature T and electron chemical potential \x e in the outer region r ~ R decrease 
to ~ 1 MeV. As a result, the disk material at r ~ R is not v-cooled anymore: R(t) 
exits the neutrino-cooled region r < r lgn on the M — r diagram (Fig. 2). The viscously 
produced heat outside r l&l is stored and advected by the spreading accretion disk. The 
spreading matter is then marginally bound to the black hole, c s ~ vk- 

o Electrons become non-degenerated. Pressure is not dominated by neutrons anymore: 
it is dominated by radiation and e ± pairs, P ~ Py + P± ~ a r T 4 . 

o Y e freezes. 

o Nuclear burning occurs: free nucleons n and p recombine into a particles. This 
process releases energy of 7 MeV per nucleon, comparable to the binding energy 
GMnipj r, and unbinds most of the disk matter, ejecting it in a freely expanding windo 

All these changes happen as R{t) grows from ~ 50 to ~ 100r g . 



A one-zone model of the spreading disk is calculated in a recent work [38]. Let us 
estimate here one characteristic radius R* at which t^ sc = ? wea k- Here f wea k is the time 
of conversion n <->• p through reactions (OD) [23I]. At i?*, pressure is already becoming 



3 The mass of the debris disk is sensitive to the parameters of the binary system before the merger, in 
particular to the mass ratio and the spins of the two companions, see e.g. 03711 for a review. 

4 In addition to nuclear burning and viscous heating, the matter is heated by neutrinos emitted at r ~ r g 
(the mass accretion rate by the black hole is still significant when R(t) approaches 10 8 cm). The energy 
deposited by neutrinos in the advective zone of the spreading disk, r ^ 50r g is comparable to the viscously 
dissipated energy in this zone. 



dominated by radiation and non-degenerate e ± pairs. On the other hand, n and p have 
not yet recombined. The timescales t wea ^ and £ v i sc are then given by 



kT \ 5 1 fH x 2 



Using the hydrostatic balance P/p — (H/r) with p « m/Aizr H, one finds 

fcr MeV, ^ « 1.0y e /i- 7 / 4 M- 3/4 mi/ 4 , (9) 

where h = 2H/R,M 3 =M/3 M , m 32 = m/10 32 g = (m/O.O5M ), and R$ =i?/10 8 cm. 
(Note that \l e jkT depends on m and h only, not R.) This gives 

^ « 17 /. 3 / 4 ao. 1 M 3 - 3/ X /2 m 32 5/4 (10) 

^visc 

=^ /?*«4.5x 10 7 /i" 3/14 a ai 2/7 M 3 3/14 m3/ 14 cm, (11) 
where GCo.i = a/0. 1 . The density and temperature of the disk at r = i? = R* are 

W; » l3hKf 3 *c$Ml /2 *n£' 2 * MeV, (12) 

p^2x l^fc-VWc^^/w^i/M g cm _ 3 _ (13) 

At this temperature and density matter is close to the neutronization line kT n (p) = 
1.04 (p/10 8 g cm- 3 ) 1 / 2 MeV (Fig. 2), 

^»0. 9 »''\W 14 % 1/H (14) 

Hence the equilibrium value of Y e at i?* is y* ~ 0.5. The actual Y e in the spreading 
disk gradually freezes out as R(t) passes through R* and its asymptotic value after 
the transition can differ from F*. The freeze-out Y e ~ 0.3 is found in |[38|l using a 
dynamical model for the spreading disk. Note that the model describes the average value 
of Y e . Viscous spreading is a random diffusion process, so different elements of the disk 
spend different times near R+, and a longer residence time at R* gives a higher Y e . One 
can therefore expect a mixture of different Y e in the spreading disk, with a dispersion 
AY e /Y e ~ 1 around the average value. 

Soon after passing R* this mixture is heated by nuclear recombination and ejected in 
a wind of duration 

r °~ J ^visc (-^ — 10 8 cm). Subsequent nucleosynthesis in the expanding 
ejecta produces diverse radio-active elements, including some with a long life-time. 
Their decay can make the ejecta visible to a distant observer Q39Q. In particular, material 
with Y e m 0.5 will synthesize 56 Ni. 56 Ni decays when the ejecta expand so much that 
their thermal radiation can diffuse out and escape to observer, producing an optical flash 
similar to normal supernovae. 



Spreading of collapsar disks 



The collapsar disks are continually fed by the infalling stellar matter during a long 
time ?i n f a ii ~ 10 s (and longer, with a decreasing infall rate). The model posits that the 
angular momentum of the infall, 4nfall> is sufficiently large to form a viscous disk 
e.g., the circularization radius of the infall in the numerical model of [|9|] is r c i rc < 30r g . 
The accretion timescale at this radius, t v - lsc ~3x 10 _2 a^s is much shorter than ?i n f a ii. 
This led |@] and many subsequent works to picture a low-mass disk, m ~ W KC M ~ 
3 x 10~ 3 M that is continually drained into the black hole and re-filled with fresh infalling 
matter. 

The picture of a low-mass viscous disk is, however, implausible. Conservation of 
angular momentum requires the following: (1) The disk spreads during ?i n f a n to a radius 
R ~ 3 x 10 2 r g (oc/0.1) _1 where the viscous timescale is comparable to ?i n f a ii - (2) The disk 
accumulates mass m that carries angular momentum 7 = J m — 7 acc . Here 7 to t = Mice Wall 
is the total angular momentum processed by the collapsar disk, 7 acc ~ M acc r g c is the 
angular momentum accreted by the black hole, and M acc ~ M^aii ~M Q is the mass 
accreted through the disk. Since 7 to t > 7 acc for any viscous disk and usually 7 to t 3> 7 acc , 
such disks must store 7 ~ 7 to t (unless almost all angular momentum is carried away by 
a wind). This implies that the disk accumulates the mass, 

m^M a J-f^-=M acc (— ) ' ~M . (15) 

'infall \ ^circ / 

The disk mass may be much smaller than this estimate only if /i n f a ii is so small that 
7tot ~ 7 acc , which leaves 7 <C 7 to t for the disk. This condition leads, however, to the 
inviscid mini-disk described in the previous section. A low-mass viscous disk could 
form only if /i n f a n is fine-tuned toward the boundary between the viscous and mini-disk 
accretion regimes (cf. Fig. 7)0 

Spreading of viscous disks in collapsars through a radius ~ 10 2 r g is accompanied 
by significant changes, similar to the evolution of merger disks described above. In 
particular, matter acquires a positive Bernoulli constant as a result of viscous, nuclear, 
and neutrino heating. The infalling material of the progenitor star exerts an external ram 
pressure on the disk, and can confine the disk initially, but eventually the disk pressure 
must win and its matter will expand with a velocity ~2x 10 9 cm/s. The ejected mass 
~ 10 33 g carries the energy E ~ 10 52 erg and will explode the outer parts of the star. 
To a first approximation, the expansion of the outer disk may be described as a thermal 
explosion driven mainly by nuclear burning of n and p into a particles. A fraction of 
the unbound disk matter will turn into 56 Ni and should create a bright supernovae-like 
event. 



The very small m ~ 0.003 M G found in the simulations of [ 9] may be the result of the imposed absorbing 
boundary condition at r m = 50 km sa 5r g , which is ~ 5 times larger than the true inner radius of the disk, 
r ms ~ r x . The large r m implies an artificially large 7 acc , which happens to be nearly equal to y tot in the 
model, permitting J <C /tot- In addition, a (small) fraction of 7 tot is carried away by the wind. 



Self -similar spreading at late stages 



The disk matter that has spread beyond ~ 10 cm is largely unbound and ejected, 
however, some matter remains bound and rotating in a remnant disk. Its mass is hard 
to estimate; it could be as large as ~ 0. 1 M Q for collapsars and ~ 0.01 M for mergers. 
This remnant is composed of recombined nucleons, and further fusion reactions are not a 
significant source of energy (compared with the virial/gravitational energy). The central 
source of neutrinos switches off as M drops, so neutrino heating is also insignificant. 
This advective remnant disk will continue to spread viscously to larger radii, gradually 
draining its mass into the black hole and possibly losing mass to a wind. 

If the mass loss through a wind is small, the spreading enters a simple self-similar 
regime such that J = const, R[t) grows as a power-law with time, while M(t) and m[t) 
decrease as power-laws with time. Detailed self-similar models of this type were studied, 



see e.g. H40L 14111 . The advective disk has a scale-height H ~ r, sound speed c s ~ vk = 
(GM/r) 1//2 , and Q. ~ Qk- Its kinematic viscosity coefficient is v ~ Ctc s H ~ av^r = afc, 
where /k(^) = (GMr) 1 / 2 . The disk spreading is a diffusion process described by 



/ t \ 2/3 

R 2 (t)~v(R)t~a(GMR) l / 2 t R(t) ~ a 2/3 (GM) 1/3 t 2/3 ~R l-j , (16) 

where subscipt "0" refers to an initial reference moment of time to. The disk mass m(t) 
is then found from the condition J = (GMR) 1//2 m = const, 

J fR\- l/2 ft\- l/3 

W (GMR) 1 / 2 l) \R J \toJ 



and the accretion rate is given by 



-4/3 



M(t) ~ - ~M ( L ] . (18) 



This self-similar solution may not apply if the disk loses mass through a wind and [38] 
consider solutions that include the wind. In general, advective disks are only marginally 
bound by the gravitational field of the black hole and their Bernoulli constant can be 
positive i42l. 143I1 . This is expected to cause a strong wind. On the other hand, bound 
solutions with a negative Bernoulli constant were found for spreading advective disks 



114 111 . The mass loss through a wind then depends on the poorly understood vertical 
distribution of viscous heating inside the disk and the behavior of the magnetic field 
above the disk. 



CONCLUSIONS 

Hyper-accretion disks are formed from matter with a modest angular momentum, with 
circularization radius r c i rc well inside 10 2 r g w 10 8 cm. Matter with r c j rc Z I0r g is sup- 
ported by the centrifugal barrier, and its accretion is driven by viscous stresses on a 



timescale to ~ 0.1(oc/0.1) 1 s, where a ~ 0.01 —0.1 is the viscosity parameter. These 
viscous disks are dense and hot, and emit copious neutrinos as long as the accretion rate 
M exceeds M ign ~ 0.03 ( a/0. 1 ) 5 / 3 M s~ 1 . 

Neutrino annihilation above the disk deposits a significant energy that can power GRB 
explosions. In contrast to previous expectations, the rate of energy deposition E V y is 
found to be not sensitive to the details of neutrino transport and the vertical structure of 
the accretion disk. It is given by E V y m Eq(M/M & s _1 ) 9//4 in a broad range of accretion 
rates M- lgn < M < Mtrap (eq.©. The normalization factor Eq is very sensitive to the black 
hole spin a. For instance, Eq « 10 52 erg s _1 is found for a black hole with a = 0.95, 
which is two orders of magnitude larger than Eq for the Schwarzschild case a = 0. 

Remarkably, all neutrino-cooled viscous disks regulate themselves to a characteristic 
state such that electrons are mildly degenerate, Y e ~ 0.1, and free neutrons dominate 
the pressure in the disk II 1711 . The neutron-rich matter may contaminate the jet from 
the accreting black hole and get ejected with a high Lorentz factor. Then the gradual 
decay of the ejected neutrons affects the global picture of GRB explosion on scales up 
to 10 cm, where the GRB blast wave is observed [|44l . 



The disk size R grows with time as a result of viscous spreading. Most of the disk 
mass m(t) resides near R(t) and its state is not described by the steady model (which 
remains valid at radii r < R). Instead, it is described by a markedly different spreading 
solution. In particular, as the disk spreads to ~ 100r^ it is heated both viscously and 
by the nuclear burning of free nucleons into helium^] As a result, the disk is disrupted 
before it spreads much beyond 10 8 cm: the heated flow acquires a positive Bernoulli 
constant and gets unbound. 

Then most of the disk mass m is ejected with a velocity ~ 0.1 c and a total energy 
~ 10 52 (m/M Q ) erg. A fraction of the ejected matter acquires Y e « 0.5, which favors the 
synthesis of 56 Ni as the ejecta expand and their temperature drops. The ensuing gradual 
decay of 56 Ni should produce a visible optical flash on a week timescale — a supernova- 
like event. The flash is expected to be especially bright for collapsars that develop 
massive spreading accretion disks. A similar (but weaker) flash should be produced by 
the spreading disks around a merged binary; the ejected mass can be a few orders of 
magnitude smaller in this case. 

56 Ni-rich matter may also be ejected from the inner, geometrical ly thin, neutrino- 



cooled disk. This can occur if the inner disk produces a strong wind Il47l . I48I . I49L 15011 . Such 
winds are modeled as quasi-steady magnetized outflows, illuminated by neutrinos which 
can heat and de-neutronize the wind material. The details of this plausible mechanism 
are uncertain because the mass outflow rate and the asymptotic Y e in the wind is hard to 
predict with confidence — both depend on the assumed MHD behavior of the disk and 
its corona. 



6 In contrast, when matter accretes through ~ 10 2 r ? (as in the steady-state model), this process is 
reversed: helium is disintegrated, which leads to cooling. 



Following the main burst, the accretion rate is determined by the amount of matter 
that remains bound and rotating around the black hole0 M(t) decreases steeply when 
the disk spreads beyond 10 8 cm and most of its matter is ejected, however, some 
matter remains bound and continues to accrete. The evolution of M may be related 
to the observed puzzling features in the afterglow emission of GRBs. The afterglow 
is likely to be produced by the relativistic blast wave driven by the jet from the central 
engine. Its luminosity is determined by the energy and magnetization of the jet as well 
as the density profile of the ambient medium at r ~ 10 15 — 10 17 cm. A long-lived jet 
of luminosity Lj et = £j et A/c 2 would certainly impact the afterglow emission. However, 
current theories are unable to reliably predict the evolution of M and £j et . For instance, 
one could speculate that the jet switches off abruptly as M decreases below a threshold, 
which causes the observed steep decay in the afterglow light curves. 

While the mechanism of the relativistic jet and its evolution with M(t) remain un- 
certain, the non-relativistic massive ejecta with v ~ 0.1c is a robust consequence of 
viscous-disk accretion. Viscous disks certainly form in merger events. The standard col- 
lapsar model also assumes the formation of a viscous disk, but this case is less certain. 
The minimum angular momentum needed to form a disk is ~ r g c, and collapsars were 
proposed as rare events of stellar collapse with I > r g c. Hence, statistically, the accretion 
flows in collapsars are likely to have small / and their disks can be smaller than ~ lOr^,. 
Such mini-disks are not centrifugally supported and accrete faster than viscous disks. In 
contrast to the viscous regime, this low-angular momentum accretion leaves no remnant 
disk in the end of the core collapse, involves no viscous spreading, and may not eject 
much mass. However, it still can produce a powerful relativistic jet via the Blandford- 
Znajek mechanism and/or neutrino heating near the rotation axis. 
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